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Abstract 

We further analyze a class of recently studied metastable string vacua obtained by 
wrapping D5-branes and anti-D5-branes over rigid homologous S'^'s of a non-compact 
Calabi-Yau threefold. The large dual description is characterized by a potential 
for the glueball fields which is determined by an auxiliary matrix model. The higher 
order corrections to this potential produce a suprisingly rich phase structure. In 
particular, at sufficiently large 't Hooft coupling the metastable vacua present at 
weak coupling cease to exist. This instability can already be seen by an open string 
two loop contribution to the glueball potential. The glueball potential also lifts some 
of the degeneracy in the vacua characterized by the phases of the glueball fields. This 
generates an exactly computable non-vanishing axion potential at large A^. 
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1 Introduction 



At present, few exact results are available on the phase structure of supersymmetry 
breaking backgrounds in either gauge or string theories. Indeed, whereas the re- 
quirement of holomorphicity determines the form of many corrections to the vacua 
of a supersymmetric theory, generically no such constraint is available when super- 
symmetry is broken. It is nevertheless natural to consider non-supersymmetric 
metastable vacua of a supersymmetric theory because the underlying supersymme- 
try of the theory allows more control over the dynamics of supersymmetry break- 
ing. Non-supersymmetric metastable string constructions have been studied in 
[D [21 El in El El El • Recent progress in finding non-supersymmetric metastable vacua 
in supersymmetric QCD-like field theories was achieved in [8] and subsequent string 
theory realizations of this work [9l [TOl [TTl [T2] . 

On the other hand, it is by now well-established that in certain cases the large 
N supersymmetric dynamics of open strings at strong 't Hooft coupling admits a 
holographic dual description in terms of weakly coupled closed strings. Notable 
examples are the AdS/CFT correspondence [T3l|lll|T5] and geometric transitions [U 
[T6l [T7] . In this note we exploit the fact that large holography is expected to be a 
more general property of many non-supersymmetric gravitational systems to analyze 
the phase structure of a strongly coupled supersymmetry breaking background. In 
particular, we further study the large A^ dual of a configuration of branes and anti- 
branes with J\f = supersymmetry of the type recently considered in [6]. 

More precisely, we study type IIB string theory compactified on the local Calabi- 
Yau threefold given by a small resolution of the hypersurface defined by: 



where W'{x) is a polynomial of degree n and x, y,u,v G C. Our brane configuration 
consists of spacetime filling D5-branes and anti-D5-branes wrapped over homologous 
and minimal size rigid S'^'s of the internal geometry. We denote by |A^j| the number 
of branes or anti-branes wrapped over the i^^ S*^. Here Ni is understood to be 
a positive (negative) integer for D5-branes (anti-D5-branes). In the absence of 
branes and anti-branes, the resulting theory in four dimensions would have preserved 
M = 2 supersymmetry. This system is non-supersymmetric because each type of 
brane preserves a different M = 1 supersymmetry. Even so, this configuration is 
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metastable because the tension of the branes generates a potential barrier against 
the expansion of the S^'s. 

In the holographic dual theory the original branes and anti-branes wrapping 
S'^'s are replaced with flux threading topologically distinct S^'s of a new geometry 
described by a hypersurface of the form: 

y2 = ^ l,^_^^n-l + . . . ^ ^ (2) 

where the bi are normalizable complex deformations of the singular geometry. The x 
and y coordinates define a Riemann surface fibered over the coordinates u and v. As 
shown in figure ([T]), each of the new Calabi-Yau reduces to a contour encircling a 
finite length branch cut of the complex x-plane. The size of each is controlled by 
a flux induced effective potential. Supersymmetric configurations of this type have 
been studied in [18]. It was conjectured in [6] that this same geometric transition 
remains valid for non-supersymmetric configurations. In this case, the vacuum is 
determined by the full potential rather than the superpotential. Rather importantly, 
as opposed to a generic A/" = 1 theory with broken supersymmetry, at leading order 
in the form of the Kahler potential is fixed by the special geometry of the 

manifolcil]. For sufficiently low flux quanta the size of each is stabilized at a small 
value [B]. Throughout this paper we will refer to this location in moduli space as 
the semi- classical expansion point. 

But as we increase the 't Hooft coupling the S^'s will expand in size so that 
higher order corrections to the effective potential will play a more important role in 
determining the vacuum of the theory. These corrections are efficiently summarized 
by an auxiliary matrix model [T91I201I2T] which determines the classical periods of the 
closed string dual. One of the purposes of this paper is to show that such corrections 
generate an intricate phase structure which is absent in the supersymmetric case. 

When supersymmetry is broken, a metastable system of fluxes will eventually 
lower its energy by annihilating flux lines. Rather than treating an individual 
system, it is therefore more appropriate to treat the totality of all possible flux 
configurations which admit metastable vacua. To this end, we invert the question 
of finding critical points of Ves and instead ask: Given a point in moduli space, 
which brane/anti-brane configurations would stabilize the moduli at this geometry? 

^As explained in [6], this is due to the fact that at leading order in 1/iV, the fluxes spontaneously 
break Af = 2 supersymmetry. 
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Figure 1: Depiction of the geometric transition from the open string picture with 
branes wrapped over minimal size S'^'s (top) to the large N closed string dual where 
the branes and homologous S^'s have been replaced by flux threading topologically 
distinct S^'s of a new geometry (bottom). The lines with red crosses denote finite 
length branch cuts on the Riemann surface of the local geometry after the transi- 
tion. In the open string picture the area of the S"^ in the middle of the bulge is 
approximately |iy(a;)|. When branes are wrapped over these minimal size S'^'s, 
the bare tension of the branes creates a potential barrier against brane/anti-brane 
annihilation. 
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Framing the question this way, we find necessary conditions for the existence of 
mutually non-supersymmetric metastable configurations of D5- and anti-D5-branes 
wrapped over the S^'s of the original geometry. We next analyze the phase structure 
of the two cut geometry. 

Except when explicitly noted, the remainder of our results apply to the two cut 
geometry with the corresponding S'^'s supported by purely RR three form fiux. By 
determining the vacua of some special fiux configurations, we argue that A^i, and 
^YM respectively control the sizes and relative orientation of the cuts. For A''i = —N2 
and ^YM arbitrary, we find metastable vacua such that the branch cuts are mirror 
refiections across an axis of the complex x-plane. For A'^i and arbitrary and 
9ym = 0, we find metastable vacua such that the branch cuts align along a common 
axis of the complex x-plane. 

Even so, both supersymmetric and non-supersymmetric fiux configurations admit 
many other potentially metastable vacua. Indeed, even non-supersymmetric U (N) 
Yang-Mills theory admits 0{N) metastable vacua which become exactly stable in the 
N = 00 limit [221 ES] . Heuristically, these vacua are non-supersymmetric analogues 
of the Tr(— 1)'^ = energetically degenerate supersymmetric confining vacua present 
in pure A/" = 1 SU{N) super Yang-Mills theory [23]. For multi-cut geometries 
these other vacua are obtained at leading order in the closed string dual by rotating 
the branch cuts by a discrete angle ir/ \Ni\. Although this geometric symmetry is 
deformed by higher order corrections, the corresponding supersymmetric confining 
vacua remain exactly degenerate in energy. 

By contrast, we find that for non-supersymmetric confining vacua, the two loop 
contribution to the potential lifts this degeneracy so that it is energetically favorable 
for the branch cuts to align along a common axis. This corresponds to an alignment 
of phases in the glueball fields. Physically this follows from the fact that the potential 
energy due to the Coulomb attraction between the branes and anti-branes is lowest 
for such a configuration. We find that the energy dependence of the vacuum as 
a function of ^ym agrees with the form expected in large gauge theorie^. In 
addition, we also find a large number of potentially metastable vacua in accord with 
expectations from large N arguments. Indeed, although there is a single energetically 

^Upon imbedding our non-compact geometry into a compact Calabi-Yau threefold, this reahgn- 
ment generates a potential for the axion. Although we do not develop this into a fully viable 
phenomenological model, we believe that this mechanism may be of independent interest for solv- 
ing the strong CP problem. 
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preferred confining vacuum, the rate of decay to this lowest energy configuration is 
suppressed by large N effects. Thus, once we have established the existence of 
a single metastable brane/anti-brane configuration, general arguments from large 
gauge theories imply the presence of a large number of additional potentially 
metastable vacua. 

Restricting further to the cases A''i = — and |A''i| ^ |A'^2| with the branch 
cuts aligned with the energetically preferred configuration along the real axis of the 
complex x-plane, we find that for sufficiently large 't Hooft coupling, but far before 
the cuts touch, the theory undergoes a phase transition which lifts the metastable 
vacua present at weak coupling. This is a novel phenomenon where strong coupling 
effects lead to a loss of stability in a classically metastable brane/anti-brane system. 

Once the 't Hooft coupling is large and metastability is lost, we can ask about 
the fate of the vacuum. In order to address this, it is necessary to go beyond the 
regime where a perturbative computation of the potential is valid. Since we have the 
exact potential at large A^, we can study this regime as well. Using a combination 
of numerical and analytic arguments, we find that the dynamics of the fluxes drive 
the moduli to a configuration where the branch cuts nearly touch. Close to this 
region in moduli space, a gas of nearly tensionless domain walls will typically cause 
the system to tunnel to a metastable vacuum of lower flux. When this does not 
occur, the cuts can touch and additional light magnetic states condense. In this 
case, we find that the resulting geometry is a non-Kahler manifold. 

The organization of the rest of this paper is as follows. In section [2] we establish 
notation and review the conjecture of [6] on the large N dual of spacetime filling D5- 
branes and anti-D5-branes wrapped over S'^'s of a non-compact Calabi-Yau threefold. 
In section [3] we derive necessary conditions for the existence of a metastable vacuum. 
Beginning in section H] we specialize to the two cut geometry and explain how the 
fluxes control the sizes and relative orientation of the branch cuts. We next show in 
section [5] that a two loop effect lifts the degeneracy in energy between the confining 
vacua of the theory so that it is energetically favorable for the branch cuts (i.e. the 
phases of the glueball condensates) to align along a common axis. In this same 
section we also find a large number of additional potentially metastable vacua and 
compare the 6ym dependence of the vacuum energy density with general expectations 
from large N gauge theories. In section [6] we show that for sufficiently large 't Hooft 
coupling, this two loop effect lifts the metastable vacua present near the semi-classical 
expansion point. This causes the branch cuts to expand until they nearly touch. 
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Section [7] discusses the behavior of the system near this region of moduh space, and 
section [8] presents our conclusions and possible avenues of further investigation. 

2 Geometrically Induced Metastability 

In this section we set our notation and discuss in more detail the large dual 
description of the metastable brane/anti-brane configuration we shall study in this 
paper. The open string description of our system consists of D5-branes and anti-D5- 
branes which fill Minkowski space and wrap n minimal size S'^'s of a local Calabi-Yau 
threefold defined by the hypersurface: 

y'^ = W\xf + uv (3) 

where x,y,u,v e C and W'{x) = g{x — ai) ■ ■ ■ (x — a„) is a degree n polynomial. 
Because the geometry is non-compact, the correspond to non-normalizable modes 
which determine the relative separation between the branes. The minimal size 
S^'s of the geometry are all homologous and are located at the points where W'{x) 
vanishes. Indeed, at a generic point of the complex x-plane the area of an S"^ is 
given by the relation: 

A{x) = (^\W'{x)f + \rfy^^ (4) 

where r denotes the size of the S*^ at x = a^. Because the branes and anti-branes 
preserve different supersymmetries, the corresponding system does not preserve any 
supersymmetry. While there is no topological obstruction to the branes and anti- 
branes annihilating, the system is nevertheless geometrically metastable because the 
bare tension of the branes produces a potential barrier against the expansion of the 
branes. See figure ([1]) for the local behavior of this configuration. 

In the holographic dual description, the branes and anti-branes wrapping n ho- 
mologous S'^'s of the original geometry are replaced by fluxes threading the n topo- 
logically distinct S'^'s of the new geometry. The local Calabi-Yau threefold after the 
transition is defined by the equation: 

y2 = W'{xf + bn-lX"~' + --- + bo + UV (5) 

where the bi correspond to the n normalizable complex deformation parameters of the 
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Figure 2: Depiction of the complex x-plane corresponding to the Riemann surface 
defined by equation ([5]) with uv = 0. The compact A-cycles reduce to counterclock- 
wise contours which encircle each of the n branch cuts of the Riemann surface. The 
non-compact -B-cycles reduce to contours which extend from x = Aq on the lower 
sheet (dashed lines) to x = Aq on the upper sheet (solid lines). 

Calabi-Yau. This complex equation defines a two-sheeted Riemann surface fibered 
over the u and v coordinates. The bi split the double roots of W^x)"^, creating n 
finite length branch cuts on the complex x-plane of the Riemann surface. See figure 
([2]) for a depiction of this geometry. 

The n 5'^'s correspond to n 3-cycles Ai such that Ai (1 Aj = for all Dual 
to each A-cycle is a non-compact S-cycle such that Ai fl Bj = —Bj (1 Ai = 6ij 
and Bi Ci Bj = for all At the level of the Riemann surface, the Ai reduce 

to n distinct counter-clockwise contours encircling each of the n branch cuts of the 
Riemann surface and the Bi reduce to contours which extend from the point x = Aq 
on the lower sheet to the point x = Aq on the upper sheet. The IR cutoff defined 
by Aq in the geometry is identified with a UV cutoff in the open string description. 
The periods of the holomorphic three form Q along the cycles Ai and Bi define a 
basis of special coordinates for the complex structure moduli space: 



where denotes the genus zero prepotential. In the absence of fluxes, each Si 




(6) 
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corresponds to the scalar component of a ?7(1) A/" = 2 vector multiplet. Once 
branes are introduced, each 5*^ is identified in the open string description with the 
size of a gaugino condensate. Defining the period matrix: 

dS, dS,dS/ ^ ' 

to leading order in the expansion, the Kahler metric for the effective field theory 
is ImTij. For future use we also introduce the Yukawa couplings: 

^ ds,ds,ds,- 

We now describe the large dual description of brane configurations which 
preserve Af = 1 supersymmetry. To this end, recall that for supersymmetric fiux 
configurations the fiux- induced superpotential is [25] : 

Weff = [ H;An = a{Si + --- + Sn) + (iViHi + ■ ■ ■ + Nnlln) (9) 



where denotes the net three form fiux after the system undergoes a geometric 
transition: 

H3 = -f^RR + TiiB-f^NS- (10) 

In the above equation, i^RR is the net RR three form field strength, H^s is the 
net NS three form field strength and tub = C*o + ^^"''^ is the type IIB axio-dilaton. 
Explicitly, 

N, = Jhs, a = a, = -jH, (11) 

for all i. Note in particular that when Ni = pi + ruBQ'i is a general complex number, 
this describes the geometric transition of a {pi, qi) 5-brane wrapped over the i^^ 
S"^ with Pi units of D5-brane charge and qi units of NS5-brane charge. In the open 
string theory, the parameter a corresponds to minus the complexified gauge coupling 
evaluated at the UV cutoff: 

« = a(Ao) = — —. (12) 

27r 9ym 
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The purely Si sector of the theory is described by the Lagrangian density: 



Cs = ^Aut (Im r),j d^S'd^^S^ + A^v^ff {S^,SJ) (13) 
9s 

where Auv is a UV mass scale which is potentially different from Aq, and Vcs is given 
bjH: 

a?nW= («fc + A^'"rfe,fcU (a^ + r^'iv''). (14) 

When it will not cause any confusion, we will work in units where Auv is normalized 
to unity. 

Having reviewed the large dual description for M = \ brane configurations, we 
now describe the A/" = analogue of this description when some of the branes are 
replaced by anti-branes. In [B] it was conjectured that the form of the flux-induced 
effective potential for the Si's, is essentially unchanged from the supersymmetric case. 
To properly compare the energy of both branes and anti-branes simultaneously, it is 
appropriate to shift V^g by a multiple of the bare tensions of the branes [6]: 

Kff ^ Feff + (A^i + ■ ■ ■ + Ar„) . (15) 
5'ym 

2.1 Matrix Models and V;ff 

In this subsection we review the connection between matrix models and special ge- 
ometry. As originally proposed in |T9j, the genus zero prepotential of the geometry 
defined by equation ([5]) is exactly computed by the planar limit of a large M auxiliary 
matrix model with partition function: 



1 



Vol([/(M)) 



y^^exp |^-^Triy($)^ (16) 



where $ is a holomorphic M x M matrix and the above matrix integral should be 
understood as a contour integral. The prepotential of the n-cut geometry near the 



^In string frame, the overall scaling of the superpotential and Kahlcr potential is Wcft/Ss and 
K/g1, respectively. When it will not cause confusion, we shall suppress the gs scaling of these two 
quantities in our computations. 
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semi-classical expansion point is given by expanding the eigenvalues of $ about the 
n critical points of the polynomial W . The usual eigenvalue repulsion term of the 
matrix model causes these eigenvalues to fill the n cuts of the geometry after the 
geometric transition. With Mj eigenvalues sitting at the z*^ cut of the geometry, 
this matrix model may be recast as an n-matrix model of the form: 



J]Vol(f/(M,)) 



/ n . 

d$i ■ ■ ■ d$„ exp V - -TrlVi^t ($i 



^MM = ^7 / u^fi ■ ■ ■ u^f^exp I 7 LiYViyi'i) iri/Kint l^^fi, 




i=l 

(17) 

in the obvious notation. We caution the reader that the numbers Mj are unrelated 
to the wrapping numbers of the branes. 

The connection between the above matrix model and the special geometry of the 
Calabi-Yau threefold defined by equation ([5]) is obtained as follows. The periods of 
the A-cycles are given by the partial 't Hooft couplings of the matrix model |19j : 

Si = gsMi. (18) 

Evaluating Zmm in the saddle point approximation, the planar limit of the free energy 
for the matrix model is identified with the genus zero prepotential for the complex 
structure moduli space of the Calabi-Yau threefold: 

] *^ measure ~l~ pert (-i^) 

where J-measure correspouds to contributions from the Vol(?7(Mj)) factors in the 
path-integral measure [26] and T^crt corresponds to perturbative contributions from 
planar Feynman diagrams: 



27rzJ'-m™,sure = X^TJ'S'f log (20) 
A An 



measure 

i=l 

n 



27r2^pert = -5^^.W^(ai) + C^.■■■^..S\'■■■S\-. (21) 



0<ii 



From the perspective of the open string theory, the contribution J^measure repro- 
duces the expected Veneziano-Yankielowicz terms in the superpotential. At lead- 
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Figure 3: Topologically distinct one loop planar diagram contributions to the pre- 
potential for the two cut matrix model. Solid lines denote branes and dashed lines 
denote anti-branes. 

ing order in the expansion of the periods about small Si, these contributions serve 
to stabilize the magnitude of the glueball fields at the exponentially small value 
~ exp (— STT^/^fyM addition to this leading order behavior, the power series 

in the 5'i's given by J-'pert produces subleading corrections to the form of the glue- 
ball potential. Although the form of such corrections are difficult to calculate for 
a general confining gauge theory, in the present case the integrable structure of the 
matrix model ensures that the Cj^...j„ are in principle calculable. 

2.2 Leading Order Behavior 

Following [6], we now review the leading order behavior of metastable critical points 
of Ves near the semi-classical expansion point. Expanding the genus zero prepotential 
to quadratic order in the S'j's, the period matrix is: 



for all i ^ j- In the above expression Aij = ai — aj is the relative separation between 
the minimal cycles over which branes or anti-branes wrap. From the perspective of 
the auxiliary matrix model, this leading order behavior corresponds to the sum of the 
measure factor terms and all one loop planar diagrams. These latter perturbative 
contributions give rise to terms in the prepotential proportional to SiSj. See figure 
([3]) for the one loop contributions to the prepotential of the two cut geometry. For 
sufficiently small \Ni\ the complex structure moduli are stabilized at exponentially 




(22) 
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Veff 




0.001 0.002 0.003 0.004gA^ 



Figure 4: Plot of Kfj/ |a| in the one loop approximation along the locus Si/ gA^ = 
— > for a flux configuration with A^i = — In this plot Aq/A ~ 10^ and 
|iVi/a| ~ 0.1. In the neighborhood of the semi-classical expansion point there is a 
single critical point which is metastable. 

small values [6]: 

S.^>o = QW"ia.)All[ (^) n ( = ) exp 

where Q denotes an A''*^ root of unity. As expected, there is a mass splitting at 
leading order between the bosons and fermions which explicitly demonstrates that 
supersymmetry is broken. Finally, at leading order there is only a single critical 
point of the physical potential near the semi-classical expansion point corresponding 
to the metastable minimum. See figure (jlj) for an example of this behavior in the 
two cut geometry. 



-lir 
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Figure 5: Two loop planar diagram contributions to the genus zero prepotential of the 
two cut geometry. Solid lines denote branes and dashed lines denotes anti-branes. 

2.3 Two Loop Corrections 

We now describe two loop corrections to Vcs- Figure ([5]) depicts the collection 
of topologically distinct diagrams which contribute to the cubic term of the genus 
zero prepotential. Rather than describe the relative contribution of each of the 
twelve topologically distinct two loop planar diagrams, we merely give an example 
of the relevant combinatorics. The combinatorial factors for the two diagrams with 
purely solid lines corresponding to the disk with two holes and two disks attached 
by a tube in figure are respectively 1/6 and 1/2. The remaining two loop 
contributions generate additional terms at cubic order in the genus zero prepotential 
and are computed in |27j . 

3 Critical Points of V^ff 

In this section we derive necessary conditions that a flux configuration must satisfy 
in order for the flux induced effective potential: 

j^j (a, + Ta'N'^ (25) 

to possess a critical point. Because such configurations are not absolutely stable, the 
amount of flux through each A-cycle will decrease over a sufficiently long time scale. 
Rather than treating one particular metastable flux configuration, it is therefore more 
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appropriate to treat the totality of all flux configurations which admit metastable 
vacua. Indeed, we ask the inverse question: Given a point in moduli space, what 
configuration of fluxes would stabilize the moduli at this point? To this end, we 
solve for the fluxes as a function of the critical point. 

In what follows, we shall not require that the flux vector be an element of 
an integral lattice. Indeed, because the critical points of the effective potential are 
invariant under an overall rescaling of a and the we may approximate the iVj 
as continuous parameters once they have been rescaled to a sufficiently large value. 
Further, we shall at first allow configurations with A'^j G C 

The critical points are given by differentiating the effective potential with respect 
to the Sj for j = 1, ■ ■ ■ , n: 

= -2zaKfr = (a + (a-r) (a + riV) (26) 

Imr Imr ^ ' 

= W (djT) V (27) 
where we have introduced two n-component column vectors V and W: 

V = -^(a + rN) (28) 
Imr ^ 

W = (a + tN) (29) 
Im r 

and for compactness we have suppressed all indices determined by matrix multipli- 
cation. Solving for the n-component vectors and a in terms of V, W and r 
yields: 

2tN = V-W (30) 

2ia = tW - tV. (31) 

It is therefore sufficient to express V and W as functions of the moduli. 

Before proceeding to the solution of these equations, we now argue that for a given 
point in moduli space there are at most 2" flux configurations such that this point is a 
critical point of Vcs- In addition to the n conditions of equation ( 1271) . we obtain n — 1 
conditions from equation (!3T!) . Indeed, because the column vector a is proportional 
to the vector with I's for all entries, the n — 1 dimensional subspace orthogonal to a 
is independent of both the moduli and the fluxes. Letting Tq, ■ ■ ■ , Tn-2 denote the 
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n — 1 independent vectors which span this subspace, the dot product of the Tj with 
equation (l3Tll yields an additional n — 1 equations: 

= W'tT, - v'rTi (32) 

for i = 0, ■ ■ ■ ,n — 2. 

Combining equations (P7|) and (152]) yields a total of 2n — 1 complex equations for 
2n complex variables. Because the critical points are insensitive to the rescaling of 
the n + 1 component row vector (a, A^i, ■ ■ ■ , Nn), we conclude that up to an overall 
rescaling by a complex number, there are a finite number of fluxes which satisfy the 
required conditions. Finally, because the system consists of n — 1 linear equations 
and n quadratic equations, the number of "critical fluxes" at a given point in moduli 
space is at most 2"^. 

3.1 Non-Supersymmetric Solutions 

We now restrict our attention to flux configurations which break supersymmetry. In 
this case, the vectors V and W each have at least one non-zero component. To 
isolate the projective nature of the solutions, we introduce affine versions of V and 
W which are completely fixed by the moduli dependent matrices r and diT. Without 
loss of generality, we may take the non-zero component of V to be Vi and that of 
W to be Wn- Now define rescaled n-component vectors u and uj such that their 
complex conjugates satisfy: 

V = ViV= Vi \V) V' = Vi{u\ (33) 

W = WnOJ = Wn \0J) W = Wn {oj\ . (34) 

where to reduce notational clutter we have switched to bra and ket notation. The 
2n — 2 non-trivial components of v and uj are therefore fixed as functions of the 
moduli by the 2n — 2 equations: 

0={uj\djT\v) (35) 
0={u\r |To) (F| r \Tk) - {V\ r |To) {u\ r \Tk) (36) 
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where j = 1, ■ ■ ■ ,n and k = 1, - ■ ■ ,n — 2 with N and a given by: 



- \N) = \u) r |To) - \cJ) {V\ r |To) (37) 
2i _ 

— \a)=T \uj) {V\ T \To) -t\u) {uj\t |To) (38) 

where C is an non-zero complex constant which is undetermined by the equations. 
Because the form of these equations is somewhat similar to those obtained in the 
study of the attractor equations of Calabi-Yau black holes [2H1 [23 EDI EI] , we will 
loosely refer to the above as our attractor-like equations. 



3.1.1 Brane Types and the Real Flux Locus 

Our goal is to study the large N dual description of a system of D5-branes and anti- 
D5-branes. Note, however, that the attractor-like equations ( !37|) and ( l38|) indicate 
that every point in moduli space is a critical point of Vcs for some configuration 
of Ni G C. To impose the additional requirement that the Ni describe D5- and 
anti-D5-branes, we must further require that all of the ratios Ni/Nj be real numbers. 
This defines an n + 1 real dimensional subspace inside the 2n real dimensional moduli 
space. 

Although a complete characterization of the real flux locus is non-trivial, a partial 
description exists in the special case when the n cuts of the geometry are all aligned 
along the real axis of the complex x-plane (so that the coefficients of the polynomial 
defining the Calabi-Yau threefold are all real) and with Aq chosen so that the Tij are 
all purely imaginary. We now show that in this case the ratios Ni/Nj are all real and 
that ^YM = 0. First note that there exists a finite neighborhood around the semi- 
classical expansion point such that Tij and J^ijk are pure imaginary. Combining this 
with the leading order behavior of the S'j's described in subsection 12.21 it follows that 
there exists a neighborhood around the semi-classical expansion point such that the 
operators {u) {uj\ and \uJ) {V\ correspond to matrices with real entries. By inspection 
of equations (1371) and (1381) . this implies that Ni/Nj is real and a/Ni is imaginary for 
all j. Note that as expected, the N^ control the sizes of the cuts and the parameter 
9ym controls the relative orientation of the cuts in the geometry. 
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3.1.2 Example: Two Cut Geometry 



With notation as in section 13.11 we have: 



Pw 1 

PvPw Pv 



\To)= \ , W) 
where we have introduced: 

= V2/V, = 
= W1/W2 

with: 



di = det 



^ 111 ^ 112 
^ 112 ^ 221 



do = det 



. V2/V1 
uj) (77| = 



Pw PvPw 

1 p^ 



WI/W2 
1 



ds ± a/c?! - 4rfic?2 

2d2 

(is ± ^/ dl — Adid2 
2di 



^ 112 ^ 221 
^ 221 ^ 222 



d-x = det 



! "3 



111 112 
^ 221 222 



(39) 
(40) 



(41) 
(42) 



■ (43) 



The ± signs of equations fj^T!) and (H2l) are correlated. The requirement that (yfy^y^ > 
leads to an unambiguous assignment of brane type for each branch. Switching from 
the + to the — branch of equations fHTj) and P2|) changes all branes (anti-branes) 
into anti-branes (branes). 



4 Fluxes and Geometry 

Unless explicitly noted, in the rest of this paper we restrict our analysis to the phase 
structure of the two cut geometry with the corresponding S^^s supported by purely 
RR three form flux satisfying the condition N1/N2 < 0. In this section we explain in 
more detail how A'^i, and Oyu respectively control the sizes and relative orientation 
of the branch cuts. 

Recall that the two cut geometry is given by the defining equation: 

y'^ = W' {xf + bix + bo + uv (44) 
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where W'{x) = g{x — ai){x — 02) and 61 and 60 control the sizes and orientations of 
the branch cuts by sphtting the double roots Cj to af. When uv = 0, this is the 
defining equation for an elliptic curve. Without loss of generality, we may take the 
Oj to be real numbers such that ai > 02- When the cuts are small, we have |18j : 

Si^ ^{at -(^lY {(^t + (^1 -(^t -(^2) (45) 
-'S'2 ~ — (a^ - 02 ) {o-t + % - '^2' - «2 ) • (46) 

It thus follows that rotating the branch cuts in the complex x-plane changes the 
phases of the S'j's. Note that equations ( l45l) and ( l46l) are corrected at higher order 
by a real analytic power series in the af. When 5*1 > > S'2, the branch cuts 
of the geometry lie on the real axis of the x-plane. Finally, for future use we set 
A = ai — 02- 

It follows from the discussion in subsection 13.1.11 that the space of critical points 
which satisfy the condition N1/N2 < defines a three real dimensional subspace of 
the four real dimensional subspace locally described by the coordinates Si and S'2. 

Although an exact characterization of this subspace is beyond our reach, we 
can still provide a crude sketch by considering various special limits. At leading 
order in the expansion of the periods, the critical points of the effective potential for 
A^i > > A^2 are [6]: 



AT- 



where the Q are N^^ roots of unity for i = 1,2 and label the distinct confining vacua 
of the low energy theory. Observe that A'^i and determine the magnitudes of Si 
and S'2. The Ni therefore determine the sizes of the branch cuts in the closed string 
dual. Further, 9ym controls the relative phases of Si and S2. Indeed, as 9ym varies, 
the branch cuts rotate in opposite directions. In this section we shall assume for 
simplicity that Ci = C2 = 1- This will necessarily limit the scope of our analysis. 
We will return to this important point later on in section [5] where we will show that 
there is an energetically preferred confining vacuum corresponding to both of the 
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branch cuts aligned along the real axis of the complex x-plane. 

In the next two subsections we show that the flux conflguration A'^i = — admits 
metastable critical points on a Z2 symmetric locus in moduli space where Si = —5*2. 
In this case the branch cuts are of equal size and are mirror reflections across the 
line halfway between ai and a2- We next show that flux conflgurations given by A^^i, 

real and 9ym = admit metastable critical points with —S2, "S*!, Aq > 0. In this 
case the branch cuts are of different sizes but are both aligned along the real axis of 
the x-plane. 



4.1 Geometry of the Z2 Symmetric Locus 

We now consider the geometry of the locus Si = —S2. Setting uv = 0, the Riemann 
surface deflned by equation (14^ is invariant under the mapping a: 

X I— i> —X, y I— i> —y (49) 

provided that g, bo, ibi G M, and Oi = —02. 

We show that invariance under this Z2 symmetry implies 5*1 = — 5*2. The S'j's 
reduce to line integrals on the Riemann surface: 

5*1 = — : (fydx, 5*2 = — : (fydx (50) 



27ri / 27ri 

where by abuse of notation we let the also refer to the reduction of the A-cycles 
to two counterclockwise oriented closed loops encircling the two branch cuts of the 
geometry. Whereas the differential element ydx is by construction invariant under 
the map a, the 1-cycles Ai and A2 transform as: 

a{Ai) = -A2, a{A2) = -Ai. (51) 

It therefore follows that: 

5i + 5^=0. (52) 
On the other hand, it follows from a general residue computation that [18] : 



^1 + ^2 = -^61. (53) 
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Figure 6: Depiction of the relative orientations of the branch cuts in the two cut 
geometry for: i) Si = —S2 and ii) 5*1 = —S2. 



We therefore conclude that when Si = —S2 > 0, the branch cuts are of equal size 
and are aligned along the real axis of the x-plane. More generally, bo (resp. bi) 
predominantly controls the size (resp. relative orientation) of the branch cuts. 



4.2 Flux Configurations of the Z2 Symmetric Locus 

We now show that there exists a finite region in moduli space around the semi- 
classical expansion point such that the flux configuration A'^i = — A^2 admits critical 
points satisfying Si = —S2. It follows from the explicit expressions for the Tij 
given in appendix A that near the semi-classical expansion point and along the locus 
5*1 = —S2- 



Til 


ri2 




~T22 


-ri2 


M - M 


' 1 


1 " 


TV2 


T22 




. -ri2 


-Til _ 


2m 


_ 1 


1 



(54) 



where M = log (Aq/A^) and: 

^111 =-^222, ^112 = -^122- (55) 

Note that the above relations do not require Aq to be a real number. 

We now apply the above relations in order to simplify the attractor-like equations 
(!37j) and (l38l) . Equation (!55|) implies that along this locus, the determinants of 
equation (l43l) satisfy di = d2 and = d^ so that the discriminant d\ — 4:did2 is a real 
number. It thus follows from the formulae in appendix A that near the semi-classical 
expansion point: 

dl - Adid2 > 0. (56) 



20 



This in turn implies: 



Pv = Pu 



(57) 



Substituting equations (IMI) and into the attractor-hke equations (1571) and 
yields: 



2i 


' iVi " 




1 


C' 


. ^2 . 


= (1 -P«.) 


-1 



2i ( . . M-M.^ 
—a = I + ri2 + pn>{ru + ^2) H ^— — (1 + 

where C is a common rescaling factor. Hence, Ni = —N2 and: 
a l + pu, f , M-M\ T22-m 

Til + T-12 H — I + 



(58) 
(59) 



N, 1 - p. 



2711 



1 - P,. 



(60) 



4.3 Flux Configurations of tiie Real Locus 



When 5*1 > > 5*2, an argument similar to the one given in the previous subsection 
establishes that in a finite neighborhood of the semi-classical expansion point, both 
and pu, of equation (HTl) and ( H2l) are again real. In this case 7^ = —Tij — {M—M)/2m 
and the attractor-like equations may be written as: 



2i 
C 



' Ni ' 




(rii 


- T12) 


+ (1- 


f Pi,P«;) 


{ri2 


- T22 


. N2 . 




. 2p,„ (ri2 


- T22) 


+ (H 


h Pt,p«,) 


(rii 


- T12 



—a = [p^puj - 1) det r H — — (1 + p„) (ri2 - r22 + Pw [tu - tu)) . 

o 2m 



(61) 
(62) 



Because Tjj — tu is pure imaginary and p„ and p^ are purely real, we conclude that 
N1/N2 is purely real. Further, when M = M, the ratio a/Ni is pure imaginary. 
We therefore conclude that for ^ym = and M = M, Ves admits critical points 
corresponding to geometries with the branch cuts aligned along the real axis of the 
x-plane. 
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5 Confining Vacua and Glueball Phase Alignment 



In this section we show that two loop corrections to our metastable brane/anti- 
brane system generate an energetically preferred confining vacuum. In the closed 
string dual this preferred vacuum corresponds to a configuration where the branch 
cuts align along a common axis. We next estimate the tunneling rate for glueball 
phase re-alignment and find the rate of decay to this lowest energy configuration is 
suppressed by large N effects which can only be countered by an exponentially small 
glueball field. Viewing our construction as imbedded inside a compact Calabi-Yau, 
string theory requires that ^ym be treated as a dynamical field. We show that 
the same effect which lifts the degeneracy between the confining vacua generates a 
potential for Oyu which is consistent with general expectations from large gauge 
theories. 

5.1 Degenerate Confining Vacua 

It is well-known that pure super Yang-Mills theory with gauge group SU (N) has 
confining vacua counted by the Witten index Tr (— 1)^ = A^. Indeed, the glueball 
field of the theory attains distinct values: 

S = CA' (63) 

where C denotes an A^**^ root of unity and A is the holomorphic scale of confinement 
for the gauge group SU{N). 

We now show that at leading order in the expansion of the periods, the confining 
vacua of the n cut geometry are also energetically degenerate. To this end, note 
that equation (1221) implies that for i ^ j, Tij is constant. The claimed degeneracy 
now follows because the critical points of V^g are given by extremizing with respect 
to the variables A^Jog {Si/ (W"{ai)Al)) for all i. 

But whereas the energy of each confining vacuum in the supersymmetric case 
is zero to all orders in an expansion of the periods, higher order corrections in the 
non-supersymmetric case should lift this degeneracy. Indeed, our expectation is that 
the Coulomb attraction between branes and anti-branes will cause the branch cuts 
of the closed string dual to align in order to more efficiently annihilate fiux lines. 
We now confirm this in the case of the two cut geometry. 
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5.2 Higher Order Corrections in the Two Cut Geometry 

At leading order, the energy density of the brane/anti-brane system is [Ej: 



e(o) = 4^ (liVii + |iV2i) - - liVii m log 

TT 



9ym 



An 



(64) 



Because E^^^ does not depend on the phases of the glueball fields, the confining vacua 
are energetically degenerate. 

At higher order the effective potential takes the form: 



VeS — Kff^ + ^crt- 



(65) 



Incorporating the two loop correction to Tjj given in appendix A lifts the degeneracy 
in energy densities: 



E = 



10 lA^il \N2\ 



IT 



(66) 



where ti = Si/{gA^) and ^2 = — 'S'2/(5'A^) are given by equations P7|l and PHI) . 
respectively. The confining vacuum with the lowest energy density is given by the 
configuration with ti and t2 as close to being real positive numbers as possible. 
Without loss of generality, the geometrical significance of this result can be seen 
when A > 0. It now follows from equations (145|1 and (146|) and the remarks below 
these equations that in this case the branch cuts are nearly aligned along the line 
joining x = ai and x = a2- 



5.3 Tunnehng Rates 

In the previous subsection we showed that the two loop contribution to Vcs lifts the 
degeneracy in energy density between the many confining vacua of the theory. In 
this subsection we compute the tunneling rate for glueball phase re-alignment. In 
particular, we find that in the strict Ni = oo limit with fixed confinement scale, these 
additional confining vacua become exactly stable. For finite Ni, this leads to the 
presence of a large number of additional metastable vacua of the type generically 
present in large gauge theories [231 121] • On the other hand, treating the glue- 
ball field as an exponentially suppressed quantity, we find that the corresponding 
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tunneling rate increases. 

In general, such tunneling events correspond to the nucleation of a bubble of vac- 
uum with lower energy density inside the higher energy density vacuum. Assuming 
that in the Euclidean continuation of Minkowski space that this bubble is an 0(4) 
symmetric configuration, the thin wall approximation of the tunneling rate is 



r ~ exp (Sbounce) = exp (~ ^2^ (AF)0 ^^^"^ 

where T is the tension of the domain wall and AV^ is the change in energy density 
between the two vacua. 

The domain wall solutions separating the confining vacua of the theory are given 
by wrapping D5-branes over the A-cycle threaded by positive flux and anti-D5-branes 
over the A-cycle threaded by negative flux. We approximate the tension of such a 
domain wall using the supersymmetric analogue with Ni,N2 > 0: 

T=- |Weff(Ci^i,C2^2) - Wefr(Ci ^1,(2^2) I (68) 
9s 

where d and C' are A^*^ roots of unity and the Si^s are evaluated at a supersymmetric 
critical poinlo. For simplicity, we now restrict our analysis to tunneling events which 
only rephase 5*1. The leading order 5*1 dependence of IIi and 112 is: 

2.m, = 5,(log^-l) (69) 

27rffl2 = 5ilog^ (70) 

The tension of the domain wall solution which interpolates between different discrete 
phase choices for 5*1 is therefore: 



T = -^\N,S,iC^-C[)\ = — 

^T^9s T^9s 



A^i^ism ^ ^ 



(71) 



where (i = exp (27rzZ/ \Ni\) and = exp (27rz/'/ |A^i|). The change in energy density 



"^Thc additional factor of l/gs in the tension formula follows from the string frame normalization 
of the superpotential described in the footnote above equation p^ . 
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between the two vacua has norm: 



\AV\ 



20|iVi| liVsl \ti 



TT 



COS 



cos 



9 + 2nl' 



(72) 



where 6/ \Ni\ denotes the argument of ti. 

We now estimate the value of T^/ (A^)'^ for different glueball phase ahgnment 
tunnehng events. In particular, we show that glueball phase alignment typically 
proceeds via a single large drop in energy density rather than a sequence of tunneling 
events with smaller drops in energy. As it is irrelevant for the considerations to 
follow, we suppress all dependence on 6 in the expressions to follow. 

First consider the tunneling event from / = A'^i/2 to /' = corresponding to a 
single instanton process with the largest possible drop in energy density. The bounce 
action is proportional to: 



Sn,/2^o a TV {AVy 



J_\N^\hl\gA 
40% lA^J^ 



3|4 



(73) 



Because the scaling of the Ni and in the large limit are all comparable, note 
that the corresponding bounce action scales as A^^ so that for finite ti this tunneling 
event is highly suppressed. Note, however, that it is also natural to consider the limit 
in which ti is exponentially suppressed. In this case the tunneling rate increases. 

Next consider the tunneling process from the vacuum / = A'i/2 to /' = A'i/2 — 5 
for 6 small compared to A^i. In this case, the bounce action is proportional to: 



S 



Ni/2^Ni/2~5 



OC TV (A^)^ 



Ai 
N2 



\ti\ \gA 



3|4 



(74) 



which scales as A^'* in the large A^ limit. The tunneling rate for a small re-alignment 
in the glueball phase is therefore much smaller than that due to a large re-alignment 
in phase. 

This conclusion is further supported by considering the tunneling process from a 
vacuum I = 6 to I' = where 6 is again small compared to |A^i|. In this case we 
find that the bounce action is given by the same expression as equation fl71|) . The 
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decay rate due to rephasing ti is therefore: 



exp 



27 1 
'¥40% 



N2 



\ti\\gA ' 



(75) 



It therefore follows that as 6 increases, the corresponding tunneling rate also in- 
creases. Tuning the parameters of the theory so that ti remains fixed, note that 
when the ratio |A'^i/A'^2| ^ 1, the tunneling rate is suppressed. Conversely, when 
the ratio |A'^2/^i| ^ 1? the tunneling rate is higher and the smaller cut aligns along 
the real axis on a shorter time scale. Physically this corresponds to the fact that 
the orientation of a cut fluctuates less as the amount of flux passing through the 
corresponding A-cycle increases. 



5.4 Axion Potential 

As we have seen above, the two loop contribution to Vcs aligns the phases of the 
glueball fields. It follows from string theory that upon imbedding our model in 
a compact Calabi-Yau threefold, ^ym must be treated as a dynamical fielcj^l with 
potential given by equation (166!) evaluated at the preferred confining vacuum. The 
effective value of the ^-angle on the branes and anti-branes follows from the relation: 

S, = QA^ = 0\A,fexp{ze,/m) (76) 

where Q denotes an lA^jl**^ root of unity. 

We now illustrate the form of the axion potential in the case A'^i = —N2 = N . 
Equations (H7I) and (HSI) imply that for such flux configurations, the phase of Aq/A 
does not contribute to the phases of the S'j's. For simplicity, we further restrict to 
the case where is purely real. In this case, the effective value of the ^-angles 
for the branes and anti-branes satisfy 9i = —62 = Oym- The energy now takes the 
form: 

20N^\t\ ( /eYM + 27r/\ /eYM + 27r/' 



where t = ti = t2 and / and /' are integers. Evaluating / and /' at the preferred 



^In a compact Calabi-Yau, the relative separation A between the branes becomes a normalizable 
mode. In this section we assume that there exists a mechanism which stabilizes this value. 
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confining vacuum configuration yields a potential for the axion: 



(0) 40iV^ |t| . / f9YM + 27d\\ AON^ \t\ fOyu 

Kx (6'ym) = mm ( cos ( — cos 



TT V ^ / 

(78) 

where in the last equality we have assumed that Oyu G [— 7r,7r). This potential 
has a minimum at 6ym = 0. In the more general case where |A'^i| ^ \N2\ and both 
A and gA acquire complex phases, the minimum of will be shifted away from 
this value. Naively, the requirement that the physics remain invariant under the 
substitution 6'ym — > ^ym + Stt appears incompatible with the 9ym/N dependence of 
the final expression of equation (ITS!) . That the physics does remain invariant follows 
from the first equality of equation (|75]) . 

In fact, this is a general issue in large gauge theories. Recall that in large N 
U (N) pure Yang-Mills theory, the Lagrangian density is: 

^ = -jyTr {F,^Fn + ^^"'"'^^^ (79) 

where A denotes the 't Hooft coupling of the theory. The dependence of the vacuum 
energy density on the 6'-angle is: 



E.ac (e) = N'f ^- j (80) 

where / is a function which is well-defined in the large limit. The factor of 
N'^ arises from the number of degrees of freedom and the 6/N dependence follows 
from the requirement that the energy density must possess a well-defined large N 
limit [331 1221 l23j. In order to restore invariance under 6^6 + 27r, it is natural to 
conjecture that the full 9 dependence of Eyac is [22] : 

E...(9)^iV=min/(^) (81) 

where the integers / label distinct metastable branches of vacua. It has been shown 
that this type of vacuum structure arises both in softly broken supersymmetric QCD 
with small gaugino masses [MJ [35] as well as from D-brane constructions of large TV 
gauge theories [23]. In the limit N = oo, it is expected that if ti remains fixed and 
in particular is not exponentially suppressed, that these additional vacua become 
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exactly stable [23l [24] . This is indeed consistent with the decay rates given by 
equations ( !73l) . (17^ and (1751) . Multiple branches of vacua may also be present in 

QCD [aniETiEHiEniiio]. 

Returning to the first equality of equation (1781) . note that our expression for the 
vacuum energy density is of precisely the form conjectured in equation flHTl) . This 
provides an exactly calculable example where the expectations discussed above are 
explicitly borne out. Indeed, although there are higher order corrections to E^ac 
and therefore to V"ax, to leading order in 1/A^, all of these corrections are captured 
by the closed string dual description. 



6 Breakdown of Metastability 

When the amount of flux through each 5*^ is sufficiently low, the corresponding 
glueball field is stabilized at a small value. But as shown in section [5l two loop 
contributions to generate a preferred confining vacuum which aligns the phases 
of the glueball fields. In this section we show that for sufficiently large values of 
the 't Hooft coupling, this same two loop effect lifts the metastable vacua present 
at weak coupling. Note in particular that because the corrections to the Kahler 
potential are of order 1/iV, the holographic dual description of the brane dynamics 
in terms of the flux induced potential becomes more accurate as the amount of 
flux through each increases. 

Before proceeding with a more precise analysis, we flrst give a heuristic derivation 
of the value of the 't Hooft coupling for which we expect higher order corrections 
to V^g to lift the metastable vacua present at weak coupling. Recall from equation 
( l64l) that the leading order energy density of the brane/anti-brane system is: 



i?(o) = 4^ {\N,\ + m) - - m m log 



Ao 



A 



(82) 



The flrst term corresponds to the bare tension of the branes and the second term 
corresponds to the Coulomb attraction between the branes. 

Returning to the discussion near equation f|T5|) . it follows that when |A'^i| > |A^2| 
and iVi > > iV2, we have: 

£;{o) > (AT^ + N^) (83) 
9ym 
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with similar inequalities for different choices of relative magnitudes and signs for 
the Ni. We expect to lose metastability precisely when the Coulomb attraction 
contribution to the energy density becomes comparable to the bare tension of the 
branes. This is near the regime where E^^^ is close to saturating inequality fl83l) . 
This yields the following estimate for the breakdown of metastability: 



' log 



9ym \N: 



A, 



A 



(84) 



where in the above expression we have dropped all factors of order unity. 

Perhaps more surprisingly, this breakdown in metastability is calculable near the 
semi-classical expansion point. Indeed, for illustrative purposes we show in the 
subsection to follow that for A''! = —N2, the two loop contribution to the glueball 
potential causes V^g to develop a local maximum. Beyond this local maximum the 
potential subsequently rolls downward beyond the regime where perturbations about 
the semi-classical expansion point provide an accurate description. 

The rest of this section is organized as follows. For simplicity and with the anal- 
ysis of section [5] in mind, we restrict to flux configurations which produce metastable 
vacua with the branch cuts aligned along the real axis of the complex x-plane. In 
subsection 16.11 we illustrate in the case Ni = —N^ that the two loop contribution 
to Kff generates a local maximum for the glueball potential. In subsection 16.21 we 
study the breakdown in metastability for flux configurations with A'^i = —N2 and 
^YM = 0, and in subsection 16. 31 we perform a similar analysis for the case |A'^i| ^ |A^2| 
and ^YM = 0. In both cases we find that the value of the flux at which metastability 
breaks down is in rough agreement with the estimate given by equation 

6.1 Two Loop Corrections to V^g 

In this subsection we show that the higher order corrections to the periods alter the 
shape of the flux induced effective potential. Along the locus 5*1 = —5*2 > 0, the 
leading order behavior of the effective potential is: 
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(2M-logt)logt (85) 
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where we have introduced the parameter t = Si/gA"^ = —S2/gA^. This potential 
possesses a single extremum which is a minimum. We now show that this behavior 
is only correct for very small t. The higher order behavior of Vcs is: 



V, 



\a\ 



At:' 



eff 



If |'(2M-logt) logt 
TT (2M + 6t - log t) I - I f I ' (68Mt + 204^2 + 28t log t) 



(86) 



To determine the appearance of a local maximum in V^fj, we compare the relative t 
dependence of all terms in the above expression. The appearance of a local maximum 
is due to the term proportional to 68Mt. We may therefore approximate V^s as: 



^ ^(0) _ 68M|iVr 

TT 



t 



2M -\ogt 



(87) 
(88) 



We now show that for a suitable range of values, this new term causes V^s to develop 
a local maximum and a local minimum. 

To this end, first consider the derivative of V^l^^ with respect to t: 
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(89) 

Note that the term in parentheses is a quadratic polynomial in logt and therefore 
achieves a maximal value. By including the contribution from dV^'^^ / dt, we find that 
for suitable values of \N/a\, the equation dVcs/dt = possesses at least two solutions. 
These solutions correspond to a metastable minimum and a nearby maximum for the 
effective potential. 



6.2 Breakdown of Metastability: Ni = -N2 

When A'^i = —N2 = N and 6ym = 0, the metastable minima of V"efr correspond to two 
equal size branch cuts aligned along the real axis of the complex x-plane. To study 
the appearance of an instability in Vcs, we return to our general approach of solving 
for the fluxes as a function of moduli. The only non-trivial moduli dependence arises 
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Figure 7: Plot of Kg/ |a| in the two loop approximation along the locus Si/ gA^ = 
—S2/gA^ > for a flux configuration with A^i = — In this plot Aq/A ~ 10^ and 
|iVi/a| ~ 0.1. 



from the ratio: 

a Atti , , 

N - 

Along the locus 5*1 = —5*2 > 0, the moduli dependent function a/N may assume 
the same value multiple times. This implies that for a given set of fluxes, there are 
multiple critical points of Ves- 

We now briefly switch perspectives and view these critical points as functions 
of a/N. When these critical points approach the same point in moduli space, the 
effective potential develops a flat direction. Viewed as a function of moduli, when 
the 't Hooft parameter: 

A = g'y^ \N\ (91) 

reaches a maximal value, the system develops an instability. 

We now determine the value of t for which A attains a maximum. It follows from 
equation ( l60l) and the expressions of appendix A that: 



-logt + M + M + t(6 + 20M + 20M- lOlogt^) +0(t^logt). (92) 
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Note that to leading order in t, A ^ ~ — logt + M + M and therefore does not possess 
a minimum. The quantity is minimized when: 

1 



(93) 



where L = log IAq/A''! + log (20e~'^^^^) and W-i is the — l**^ branch of the Lambert 
W-f unction. 

It is also of interest to study the A dependence of t^. For sufficiently small 
and A, we have: 



log 



An 



20U 



(94) 



Fixing the value of t, we conclude that the system always develops an instability 
once A is sufficiently small. In particular, we see that by tuning A to a sufficiently 
small value, can be made arbitrarily small, justifying the two loop approximation. 

To conclude this subsection we estimate the value of (^ym |A^I which Ves devel- 
ops an instability. It follows from equations fl92|) and flMl) that for sufficiently small 
A and L: 
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where A^, is the value of the bare 't Hooft coupling at the phase transition. If we 
now drop the factors 20 and log 1^1 in the above expression, we obtain: 
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up to factors of order unity. This is in agreement with the analysis leading to 
equation (|8^ . Even so, equation fl96l) should be viewed as a crude upper bound. 



6.2.1 Masses and the Mode of Instability: Ni = -N2 

Although the above analysis establishes that Ves develops an instability for suffi- 
ciently large values of the flux, it does not directly indicate the mode of instability. 
It also does not address whether additional modes of instability appear before reach- 
ing such a flux conflguration. We now show that all the other modes are stable up 
to this point and that the unstable mode of the system corresponds to the cuts re- 
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maining equal in size and expanding towards each other. To estabhsh this, we now 
compute the bosonic mass spectrum for fluxes close to the value where we expect to 
lose metastability. After canonically normalizing the kinetic terms of 5*1 and 5*2 and 
expanding Vcfi to quadratic order, the masses squared are: 

mlA = -r^ + 2a |iV| f-^^ + ri ) (9^) 

1-v \l + v^ (1 - f)7rlmriiy 

/in 7 \ 

mis = + 2a \N\ + — (98) 

mjs = -—^--^ + 2a\N\ { + -±- 1 (99) 




<A = - ^ + 2a|iV| + - (100) 



where in the above, RA denotes the real anti-symmetric mode corresponding to both 
S'j's real with one cut growing while the other shrinks, RS denotes the real symmetric 
mode corresponding to both S'j's real with both cuts growing in size together, and 
IS and I A are similarly defined for the imaginary components of the S'j's. Further, 
we have introduced the parameters: 

\N\ Imri2^ 

a = ; , V = ; (101) 

2iTt Im Til Im Til Im T22 

where all moduli dependent functions are explicitly evaluated at the critical point 
determined by the 't Hooft coupling A. In equations fl97til00p . the term proportional 
to corresponds to the leading order contribution to the masses squared computed 
in [6], and the term proportional to 2a |A^| corresponds to the two loop correction 
to this value. As expected based on general symmetry arguments, we find that as 
a function of \N/a\, m^^ approaches zero as the flux approaches the value given by 
equation (p6|l . In figure ([8]) we show the behavior of and the next smallest 
mass squared function of \N/a\ near the regime where metastability is lost. 

It is also of interest to consider the difference in masses between the bosonic and 
fermionic fluctuations dictated by the underlying M = 2 structure of the theory. 
We find that the masses of the fermions naturally group into two sets of values. 
At leading order in 1/A^, the M = 2 supersymmetry of the theory is spontaneously 
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Figure 8: Plot of the two smallest masses mj^g and function of |A^/a| for 

Ao/A ~ 10^. At a value of \N/a\ ~ 0.122 the system develops an instability. This 
value is in rough agreement with the estimate of equation fl^ . 



broken. This indicates the presence of two massless goldstinos. Labeling the 
fermionic counterparts of the gauge bosons and the S'j's respectively by ip^^ and ip^^ 
we find that when A^^i = — the non-zero masses of the canonically normalized 
fermionic fields are all equal and given by the value: 

where as before, the first term corresponds to the leading order mass and the second 
term is the two loop correction to this value. In figure ([9]) we compare the masses 
squared of the bosonic and fermionic fluctuations. 

Although we do not include the details here, we find more generally that for 
vacua which satisfy 5*1 = —52, the system develops an instability at a similar value 
of \N/a\. In this case, the mode of instability causes the cuts to expand in size and 
rotate towards the real axis of the complex x-plane. This is in agreement with the 
physical expectation that the flux lines annihilate most efficiently when the branch 
cuts are aligned along the real axis of the complex x-plane. 
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Figure 9: Plot of log(mV(|a|^ A^y)) as a function of |iV/a| for Aq/A ~ 10^. The 
bosonic and fermionic masses squared are given by equations fl97til00l) and (11021) . 
respectively. 



6.3 Breakdown of Metastability: |iVi| > |iV2| 

We now study the behavior of Vea for flux configurations with \Ni\ ^ |A'^2| and 
^YM = 0. More precisely, we also take A'^i small enough so that the two loop 
approximation of is valid. In this case, the modulus ti = Si/gA^ fluctuates 
much less than t2 = — S'2/5'A'^. Further, because the behavior of Ves is relatively 
insensitive to the value of a/N2, it is sufficient to fix a/N2 and determine the value 
of a/Ni for which Ves develops a fiat direction. 

To this end, we employ a strategy similar to that of subsection 16.21 and use the 
attractor-like equations to treat a/Ni as a function of the single modulus t2- Because 
a larger amount of fiux passes through the corresponding to ti, it is enough to 
expand a/Ni to first order in ti and to leading order in logt2- Extremizing a/Ni 
with respect to t2, we find that Vcs develops a fiat direction when: 

Svr^ M 

^ 80^' ^ ' 

Using equation (1471) to approximate the value of ti yields the value of Ai = gyu 
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for which the system develops an instabihty: 



Stt- 



1 logi^r ^ 

80 |Ao|2 
I A I 



(104) 



where Ai,* denotes the bare 't Hooft couphng at the phase transition. To obtain a 
crude estimate of when we expect to lose metastability, we treat the left hand side 
of equation fll03p as an order one number, obtaining: 
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This estimate is in accord with equation ([5 

As in subsection 16. 2[ it is important to compute the masses squared of the bosonic 
fluctuations at the metastable minimum in order to determine the mode of instability 
for this flux configuration. In this case, the appropriate linear combination of fields 
which diagonalizes the mass matrix is somewhat messier and we defer the details of 
this computation to appendix B. We find that the unstable mode corresponds to the 
smaller branch cut increasing in size at a much faster rate than its larger counterpart. 



7 Endpoints of a Phase Transition 

For sufficiently large 't Hooft coupling, the metastable vacua present at weak coupling 
cease to exist. The moduli subsequently roll to larger values so that a perturbative 
expansion in the glueball fields Si is no longer valid. As the branch cuts increase in 
size and meet each other, the 3-cycle Bi — B2 reduces to zero size. We have checked 
numerically that over this range the potential attains a minimal value only once the 
cuts touch. Because no flux passes through B1 — B2, the moduli will not be stabilized 
away from the corresponding conifold point. Near this region, the contribution of 
new light states associated to a D3-brane wrapping Bi — B2 and an instanton gas 
of D5-branes all wrapping the same cycle will determine the low energy dynamics of 
the theory. 

Let us first consider the contribution due to the D5-branes. From the perspective 
of the spacetime, a D5-brane wrapping B1 — B2 corresponds to a domain wall solution 
which separates vacua with distinct values of the flux. Indeed, for vacua near the 
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semi-classical expansion point, this is the primary mechanism by which the vacuum 
can decay to a supersymmetric flux configuration [6]. As Bi — B2 collapses, the 
corresponding action for quantum tunneling tends to zero and a gas of D5-branes 
will scan over all available flux conflgurations. This will necessarily change the shape 
of the potential for the moduli. When enough flux has been annihilated, it is then 
possible for the moduli to subsequently either quantum mechanically tunnel through 
moduli space or (if the shape of the potential has changed enough) classically roll 
back out to the semi-classical regime. Letting the vector N^, denote the critical value 
of the fluxes for which metastability is lost, the flux vector for the new metastable 
minimum will differ from A'^,, by a flnite amount. 

Next consider the contribution due to the D3-brane. As the 3-cycle Bi — B2 
collapses, the geometry approaches a conifold point such that both the Kahler metric 
and Vcff develop singularities. Just as in the supersymmetric case, the additional 
light states which resolve this singularity correspond to a D3-brane wrapped over 
Bi — i?2- Because an analysis of a generic flux conflguration would require a fairly 
precise knowledge of the periods near this region in moduli space, we restrict our 
discussion to the case A^i = —N2 with the branch cuts aligned along the real axis of 
the complex x-plane. In subsection 17.41 we present evidence that if we ignore the 
presence of nearly tensionless domain walls, the light magnetic states of the D3-brane 
condense and cause such a flux conflguration to transition to a non-Kahler manifold. 

The rather different physical nature of the D5-brane instanton gas and massless 
states contributed by the D3-brane indicate that a more detailed analysis of the end- 
point of this phase transition is likely to be quite difficult. Our general expectation, 
however, is that the contribution due to the D5-branes will typically cause the moduli 
to relax back to the semi-classical expansion point. Indeed, in the case |A^i| ^ \N2\, 
the N2 units of flux can be treated as a probe of the background flux conflguration 
determined by A''i. In this case it is doubtful that such a small perturbation could 
cause the geometry to undergo a transition to a non-Kahler manifold. 

The situation is less clear when Ni ~ — As we show in subsection 17.11 when 
A''i = —N2 = N and along the locus 5*1 = —S2 > 0, the value of the potential when 
the cuts touch is independent of A^. This implies that while the gas of tensionless 
domain walls can still drive the system back to a metastable vacuum, it is not 
energetically favorable to eliminate a small amount of flux. The additional light 
magnetic states due to the D3-brane can then potentially influence the endpoint of 
the phase transition. 
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The rest of this section presents additional details concerning the flux configura- 
tion A'^i = — = N and is organized as follows. In subsection 17.11 we show that 
when 5*1 = —5*2 > 0, the value of Vcs at the point where the cuts touch is indepen- 
dent of A^. This computation allows us to determine in subsection 17.21 the minimal 
drop in flux necessary to tunnel back out to the semi-classical regime. In subsection 
17.31 we describe in more detail the singular behavior of Ves by passing to a new basis 
of special coordinates. Using this dual magnetic description we show in subsection 
17.41 that in the absence of D5-brane effects, a condensate of light magnetic states 
causes the geometry to transition to a non-Kahler manifold. 



7.1 Energy Near Cuts Touching 

We now show that as the cuts touch along the locus 5*1 = —5*2 > 0, the effective 
potential approaches a constant value independent of A^. By virtue of equation ([5] 
Vcs may be written as: 



VeS 1 M-M 

' At-2i( det T + —Tell 



2 -r, , ^ oMziM I V 27li 



N 
a 



(106) 



where r++ = rn + r22 + 2ri2 is the "center of mass" coupling, Ten = Th + r22 — 2ri2 is 
the complex structure modulus of the elliptic curve, and as before, M = log (Aq/A^). 
Note that as the cuts touch, Teu tends to zero. Because rn = r22, this implies tu 
tends to rn and thus det r tends to zero. We therefore conclude that V^eff approaches 
the value: 

4rn + 2^ 

Next, recall from appendix A that the behavior of th is of the form [B — M)/27ii 
where i? is a function of the moduli which is independent of a, N and M. The value 
of Vtouch is then: 

Kouch = (108) 



5t + log|^ 



2 

A"' 



where Bt denotes the value of B when the cuts touch. This is manifestly independent 
of A^. Finally, we note that as A tends to zero, so too does Vtouch- 
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Figure 10: Depiction of Kfr/lal^ as a function of 5 = a]*" — = — 03 > 
with A'^i = —N2 = N . The value of V^g at the point where the cuts touch is 
independent of A^. When is small (dashed red), the system possesses a metastable 
minimum with energy density V^cL < Kouch- For intermediate values of which still 
admit a metastable minimum (dot dashed blue), the energy density of the vacuum 
is Kaeta > Kouch- For large enough values of (solid purple), the system undergoes 
a phase transition and no metastable minimum exists. 
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7.2 Flux Hysteresis 



Close to the region in moduli space where the cuts touch, a gas of nearly tensionless 
D5-branes can cause the system to tunnel back out to a metastable configuration 
with lower flux. This tunneling will cause the flux to undergo hysteresis. In this 
subsection we determine the minimal size of this jump in the flux number. Turning 
the discussion around, this also determines the range of values for which a metastable 
vacuum near the semi-classical expansion point can tunnel to the region where the 
cuts nearly touch. 




^ 6 



touch 



Figure 11: Depiction of Kfj/ \af as a function of 6 = af — = — a2 > with 
A'^i = —N2 = A^. The flgure shows flux configurations which possess a metastable 
minimum (dashed red) and those which do not (solid purple). Close to the region 
in moduli space where the cycle Bi — B2 collapses, the presence of nearly tensionless 
domain walls can cause the amount of flux to jump. To reach a metastable minimum 
near the semi-classical expansion point, the moduli then either classically roll (a) or 
tunnel through moduli space (b). 
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We begin by characterizing the possible transition points for tunnehng processes 
along the locus Si = —S2 > 0. In addition to the metastable vacua near the semi- 
classical expansion point, there is another vacuum with identical energy where the 
branch cuts overlap almost completely. To increase the scope of our discussion, we 
briefly consider more general configurations such that > > 5*2- With notation 
as in section HI the endpoints of the branch cuts in the semi-classical regime satisfy 
af > tti > a2 > a2 with special coordinates defined by the integrals: 



"1 "2 

Si = ^[ydx S2 = -^[ydx. (109) 
2m / 2m 



where the choice of signs is dictated by branch cut considerations. In the region 
where the branch cuts overlap, we instead have > > aj~ > with special 
coordinates defined by the integrals: 

S[ = ^ [ydx S'2 = --^ [ydx. (110) 
iTTz J 2m J 



It therefore follows from geometric considerations that there is a highly non-trivial 
duality in the low energy effective field theory between vacua with Si small and vacua 
with S'i small. In addition to these two physically indistinguishable configurations, 
there is the local singular minimum where the cuts touch. 

The system can only tunnel from the region where the cuts touch to a metastable 
vacuum with small Si when: 

Kncta < Kouch (HI) 



where Kneta denotes the value of V^g evaluated at such a metastable critical point 
and Vtouch is given by equation (11081) . Approximating Kneta by equation the 
flux must therefore satisfy the bound: 



^YM |iV| > (1 - J I ■ (112) 

log 1^1 V ^B, + \og\^\ 
It now follows from equations fl96|) and (11121) that in jumping from a flux configuration 



41 



which does not admit a metastable vacuum near the semi-classical expansion point 
to one which does, the flux drops by an amount: 



6\N\> — 2- (113) 

^^^Miog|f rv^t+iog|f I' 

Turning the discussion around, the range of fluxes for which it is possible to 
tunnel from a metastable vacuum near the semi-classical expansion point to the 
region where the cuts touch is: 



-^T^ ^ ^YM N > --—^ 1 - \ (114 

where the crude upper bound follows from equation (!96l) and the requirement that a 
metastable minimum exists. Note in particular that the admissible range of values 
of the 't Hooft coupling which permit such a process is parametrically tied to the 
value for which metastability is lost. Thus, nearly as soon as we increase the 't 
Hooft coupling to a value where this quantum effect can contribute, it is overtaken 
by classical effects. 



7.3 Dual Magnetic Description 

Near the region in moduli space where the cuts touch, the flelds 5*1 and S2 have 
become large and it is appropriate to change to a dual magnetic basis of flelds. 
From the perspective of the geometry, this corresponds to performing a change of 
basis which preserves the intersection pairing of the geometry. Using this dual basis, 
we now show that the effective potential develops a cusp when the cuts toucli^. 

The appropriate change of basis is dictated by the geometry of the Riemann 
surface. Along the entire locus considered, 5*1 + 5*2 has remained zero. Further, 
the cycle Bi — B2 is close to zero size. This implies that the new A-cycles are 
Ai = B2 — Bi and A2 = Ai + A2. Dual to these are new 5-cycles Bi = Ai and 
B2 = B2. Note that this new basis preserves the intersection pairing of the geometry. 
The dual magnetic coordinates Si and Ilj are therefore related to the original special 

^Recall that a difFcrcntiablc function f{x) of a single real variable is said to have a cusp at the 
point a if f'{x) ±00 as a; ^ a^. A similar definition holds for functions of several variables. 
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coordinates by: 



Si=U2- Hi 
III = Si 



S2 — Si -\- S2 
Ha = Ha. 



(115) 
(116) 



It follows from Picard-Lefschetz singularity theory that the only non-trivial mon- 
odromy arises from the transformation Si 1— > e'^'^^Si. We therefore conclude that the 
leading order behavior of II2 is regular and Hi depends logarithmically on Si. 



;ii7) 



For a general flux configuration, the induced superpotential in the new coordinates 
is: 

Weff = j H3AQ = aS2 + {N+ + iV_) Si + 2N+II2 (118) 

where N± = {N2 ± A^i) /2. Note that Wcfr is independent of IIi. In particular, 
when A'^i = — A^2 = N, the superpotential is independent of both of the Ilj's. The 
entries of the new period matrix are: 



dS,' 



(119) 



Because II2 is regular, we may approximate the Kahler metric as: 

2 
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where the Cij are non-zero constants. In the new coordinates, the effective potential 
is: 

/- ~ \ . . 1 r -A^ 1 ^ ^ 

121 



V,s(Si,S2) = [-N a] 



Im T 



-N 

a 



in the obvious notation. It therefore follows that Vcs attains a minimum at 5i = 0. 
Replacing the Si by fields Xj with canonically normalized kinetic termCl, the 



''Although such a field redefinition will in general introduce anomalies into the Lagrangian den- 
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derivative of Vcs with respect to Xi is: 

dX^ ^ S,fiS,,S2) 

where / contains at most log 5*1 type divergences. We therefore conclude that the 
potential has a cusp at 5*1 = 0. 

7.4 A Further Phase Transition 

The singular behavior of the effective potential implies the presence of additional 
light states which have been integrated out. As in the case of the conifold [HI B2] . 
these light degrees of freedom correspond to a D3-brane wrapping the vanishing 3- 
cycle Bi — B2. In A/" = 2 language, this corresponds to a hypermultiplet which is 
charged under the gauge boson of the U{1) vector multiplet with scalar component 
Si. The superpotential is now: 

Weff = aS2 + + iV_) 5i + 2iV+n2 + HSiQlQr (123) 

where /i is a Yukawa coupling and Ql and Qr denote the TV = 1 chiral multiplets 
of the hypermultiplet. At a critical point of V^gf, the Q's condense with expectation 
value: 

{QlQr) = J. (124) 

This condensate signals the presence of a new holomorphic 2-cycle S = dAi in the 
geometry with size \N/ h\ . Perhaps surprisingly, the resulting manifold is non-Kahler! 

To see this, let us suppose to the contrary that the new geometry is a Calabi-Yau 
threefold. In this case, the absence of any other normalizable (1, 1) forms implies 
that the (1,1) form k-^ which measures the volume of S determines the local metric 
of the new geometry. It now follows from Stokes' theorem that: 




(125) 



Ai S 



This implies that k^, is not closed. We therefore conclude that the resulting mani- 
sity, this will not alter the conclusions of our analysis. 
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(122) 



fold is not a Calabi-Yau threefold, but instead belongs to the category of generalized 
Calabi-Yau threefolds [13]. While it is doubtful that the resulting physical configu- 
ration is supersymmetric or even metastable, a proper analysis is beyond the scope 
of this paper and we defer a full study of this question to future work. 



8 Conclusions 

In this paper we have studied the phase structure of a strongly coupled supersym- 
metry breaking configuration of D5-branes and anti-D5-branes wrapped over ho- 
mologous rigid S'^'s of a non-compact Calabi-Yau threefold using the large dual 
description of this system. In much of this paper we focused on the closed string dual 
geometry with two branch cuts. Even in this simple case, higher order corrections to 
the potential for the glueball fields generate an elaborate phase structure which can 
already be seen at the two loop level. Near the semi-classical expansion point, this 
two loop effect lifts the degeneracy in energy density between the many confining 
vacua of the theory. When the scale of confinement is not exponentially suppressed, 
this generates a large number of additional metastable vacua. For sufficiently large 
values of the 't Hooft coupling this same effect also lifts the metastable vacua present 
at weak coupling. After this phase transition, the branch cuts expand in size until 
they are close to touching. Although the presence of nearly tensionless domain walls 
close to this region of moduli space will most likely cause the system to relax to a 
metastable vacuum of lower flux, the presence of new massless states when the cuts 
touch may also allow the geometry to transition to a non-Kahler manifold. We now 
discuss some implications of this work. 

As the separation A between the branes decreases, the glueball potential devel- 
ops an instability. Although it is tempting to identify this instability with the open 
string theory tachyon, there is a potentially serious problem with this interpreta- 
tion. Indeed, it follows from equation (p6|) that the effective potential develops an 
instability when: 



log 



Ar 



On the other hand, the tachyonic mode of the brane/anti-brane system is independent 
of because this constant factors out of all relevant open string amplitudes. We 
therefore conclude that an identification of the two instabilities is not naively correct. 
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We have also seen a preliminary indication that the number of critical points of 
Vcs crucially depends on the amount of flux in the closed string holographic dual. 
As this amount of flux changes, a local maximum and minimum may merge. Such a 
change in the number of metastable minima cannot be detected by a Morse-theoretic 
index. This may have implications for recent attempts to count the number of super- 
symmetry breaking vacua in flux compactifications. Indeed, many of the techniques 
developed thus far rely on similar indices to count the number of admissible vacua 
|44j . It would therefore be interesting to determine whether such methods properly 
account for the metastable vacua studied in this paper. 

Although a more detailed analysis of the phase structure near the region of moduli 
space where the branch cuts touch will most likely be difficult, we have seen that 
when A^^i ~ —N2, the geometry may undergo a further phase transition to a non- 
Kahler manifold. Even so, the decay of the vacuum due to flux line annihilation may 
obstruct this intriguing possibility from contributing to the phase structure of the 
theory. It is likely that a proper description of the effective theory near this region in 
moduli space will require a more generalized effective potential which treats both the 
moduli and the fluxes as dynamical variables. Even if the geometry can transition 
to a non-Kahler manifold, the resulting configuration is unlikely to be stable. It 
would be interesting to determine the endpoint of this further phase transition. 

In much of this paper we restricted our analysis to the two cut geometry. While 
this should provide an adequate characterization of "two body" interactions, it is 
possible that the interaction of three or more cuts could lead to further novel phases. 
Although we still expect the phases of the glueball fields to align in an energetically 
preferred configuration, the relative orientation between the cuts will depend on the 
location of the Oj. Indeed, treating the branch cuts as small dipole moments, an 
energetically preferred configuration may be frustrated for a large lattice of cuts, 
much as in the two dimensional Ising model on a triangular lattice. It is well-known 
in the setting of condensed matter systems that magnetic frustration can produce 
novel phases such as spin liquids and glasses. It is therefore likely that a similarly 
rich class of phenomena are present in metastable multi-cut geometries. 
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Appendix A: Two Cut Semi- Classical r^j and Tijk 

In this appendix we collect explicit expressions for Tij and J^ijk = diTjk for the two 
cut geometry defined by equation ([5]). Setting A = ai — 02, = Si/g/S.^ and 
^2 = — 5'2/5'A^, we have 



A2 

2'KiTn = logti - log ^ + (4ti + iota) + (32^2 + 182^1^2 + UStl) + 0(1^^) (127) 

27rzri2 = - log + (-lOti - 10^2) + (-91ti - 236tit2 - 91tl) + 0{t^) (128) 
A^ 

27rzr22 = logt2 - log-| + (4^2 + lOti) + {32tl + 182^1^2 + 118t?) + 0{t^), (129) 



and: 

27iigA^J^iu = - + 4 + (64ti + 182^2) + 0{t^) (130) 

2mgA^J^ii2 = -10 + (-182ti - 286^2) + 0{t^) (131) 
27rz^A3j^i22 = 10 + 236ti + 182^2 + 0{t'^) (132) 

27ngA^J^222 = -- - 4 - 64^2 - 182ti + 0{t'^). (133) 

Appendix B: Mass Spectrum for |A^i| ^ |A^2| 

In this appendix we compute the bosonic and fermionic masses for flux configurations 
with |A^i| > |A''2|, 6ym = and 5*1 > > 6*2. With the kinetic terms of the 
Lagrangian density canonically normalized, the 4x4 bosonic mass squared matrix 
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^Bosonic takes the block diagonal form: 



" ''Bosonic 



(134) 



where the A^^'^'^ are 2x2 matrices of the form: 



( 
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(RJ)_q(RJ) 



(135) 
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and A^^'^ (^'^^^) corresponds to the mass matrix for the real (imaginary) components 
of the Si's. In the above we have defined: 
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Im Tjj d Im Sj ' 
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(136) 



and for future use we also introduce: 

l^il 



27iti Im Til 



27rt2 Ini T22 



(137) 



In the above expressions the components of Tij correspond to their value at the critical 
point of Vcs and should therefore be treated as constants. When \Ni\ » \N2\, the 
masses squared and eigenmodes of the block A^^'^ are: 
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a^-2a{ lOlA^al - 



2|A^i| + 5|iV2| 
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(0,0), (138) 



(0,0), (139) 
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and the masses squared and eigenmodes of the block A^^^ are similarly: 
6^ 



m 



m 




1 - . 




1 + ^ 




1 + ^ 




1 - . 





1 ©(0,0)^ (140) 



9 9 / , , 2|iVi| +5|A^2l\ / /l + ^A' \ . ^ 

U-«- + 2a(lO|/V.K ^' j ^-^^.lj^®(0.0)„, (141 



Grouping the fermions according to the supermultiplet structure inherited from the 
A/" = 1 supersymmetry of the branes, the non-zero fermion masses are: 

1 / 2|iVi| +5|iV2| + 10|iVi|isi^\ 
^ 1 - \ ImriiTT / 

1 / 2|iV2| +5|iVi| + 10|iVi|i2i^\ 
1 — f \ Im r227r / 

with similar notation to that given above equation fll02p . By inspection of the above 
formulae, we see that the two loop correction increases the difference between the 
bosonic and fermionic masses already present at leading order. 

Keeping I -^2! fixed, we now determine the mode which develops an instability 
as the 't Hooft coupling gyuWA approaches the critical value where the original 
metastable vacua disappear. The determinant of each block of the mass matrix is: 



9YMtltl (logti logt2 - log^ (logti + logt2; 



^(R) ^ 40967rMog^ilogt2 

logti logt2 - 20ti (logti)' (log ^ - logt: 
-20t2 {loghf (log ^ - logt2) + ■ ■ ■ 

^^^^(7) 40967r«logti logt2 

9yaA4 (log^i log ^2 - log ^ (logti + logt2 
logti logt2 + 20ti (logti)' (log ^ - logt: 
+20t2 (logts)' (log^ - logt2) + ■ ■ ■ 



6 



6 



(144) 

(145) 
(146) 

(147) 
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It follows from the last line of each expression that only m^^g_^ or rn^^g^ can vanish. 
Furthermore, because Tn\^g^ » rn^^g^, the mode of instability will cause the smaller 
cut to expand towards the larger cut. For \Ni\ ^ |A''2| this occurs at a value of ti 
given by: 

Ao 



20 



log ti + log 



A 



(148) 



Note that this is similar in form to the critical value of the moduli found in equation 
(1911) for the case Ni = -N2. 
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